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Abstract. The Kepler mission has provided high quality light curves for more than 
2000 eclipsing binaries. Tertiary companions to these binaries can be detected if they 
transit one or both stars in the binary or if they perturb the binary enough to cause 
deviations in the observed times of the primary and secondary eclipses (in a few cases 
both effects are observed in the same eclipsing binary). From the study of eclipse timing 
variations, it is estimated that 15 to 20% of the Kepler eclipsing binaries have close-in 
tertiary companions. I will give an overview of recent results and discuss some specific 
systems of interest. 


1. Introduction 

In an isolated, detached eclipsing binary (EB), the eclipses should be strictly periodic, 
with a constant interval of time between successive eclipses. In these cases, the times 
of eclipse are described by a simple linear ephemeris: 


Tmw(E) - To + PbinE (1) 

where Pbin is the binary orbital period and E is the cycle number. The residuals derived 
from a fit to a linear ephemeris form an “Observed minus Computed” or 0-C diagram. 
If one measures eclipse times (ETs) for both primary and secondary eclipses, then both 
types of events may be put on a common system by means of a phase offset 5 that is 
applied to the cycle numbers of the secondary eclipses. A single period can then be fit 
to all of the ETs. We will refer to the resulting 0-C diagram a “Common Period 0-C” 
or CPOC diagram. 

There are a number of situations where the intervals between successive eclipses 
are not constant, thereby leading to eclipse timing variations (ETVs). The points in 
the 0-C or CPOC diagrams will no longer be scattered about the horizontal axis, and a 
model that goes beyond a simple linear ephemeris will be needed. We briefly discuss 
mechanisms for ETVs that apply to stars that are well within their respective Roche 
lobes, where mass transfer and mass loss can be neglected. 

If the EB is part of a triple system, then the eclipses will either be early or late 
owing to light travel time (ETT) changes as the EB moves about the center of mass of 
the triple system. In these cases, the ETs are no longer described by a simple linear 
ephemeris (e.g. Irwin 1952): 


Tmm{E) — Tq+ PbmE + K 
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where is the eccentricity and ojj is the longitude of periastron of the tertiary orbit, V 3 
is the true anomaly in that orbit, and where the semi-amplitude K is given by 


K - 


a 12 sin/^ 
c 


(3) 


where an is the semimajor axis of the EB orbit about the third star, is the inclination 
to the observer’s line-of-sight, and where c is the speed of light. The minimum mass of 
the third body can be computed in a way that is similar to the well-known result that 
can be derived from a radial velocity curve: 


/(Ms) - 
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where Mi and M 2 are the masses of the stars in the EB, M 3 is the tertiary mass, and 
G is the constant of gravitation. Many close EBs are found to have third components 
based on ETT studies (e.g. Tokovinin et al. 2006). 

If the third star is in a relatively close orbit with the EB, then the changing prox¬ 
imity of the third body can affect the orbital period of the inner binary. Borkovits et 
al. (2003) presented an analytic study of the ETVs caused by an exterior, perturbing 
object on a large, eccentric orbit. The signals one sees in the 0-C diagrams can be quite 
complex, depending on the eccentricities of the EB and outer orbit, and on the mutual 
inclinations of the orbits. 

Before the launch of Kepler, EBs where dynamical elfects dominate the ETVs 
were comparatively rare. One example is lU Aurigae (Ozdemir et al. 2003), which is 
a massive binary (Mi = 21 . 4 M 0 and M 2 = 14.5 M©) with an orbital period of 1.8 
days. There is a third body with a period of 293 days which causes precession of the 
inner binary on a time-scale of w 335 years. If the ETVs are taken to be from ETT 
effects alone, then the minimum mass of the third body is about 15 M©, which is much 
too massive, given the amount of third light measured from spectroscopy and from the 
light curve solutions. Hence Ozdemir et al. (2003) argue that a large part of the signal 
seen in the 0-C diagram is due to dynamical effects. Another example is SS Eac, which 
was known to be a deeply eclipsing EB around the year 1900, but then stopped eclipsing 
around the year 1950, presumably due to the influence of a third body. This outer body 
with a period of 697 days was detected spectroscopically by Torres & Stefanik (2000). 
In a subsequent paper, Torres (2001) modeled the depth changes of the eclipses and 
found a precession period of about 600 years. More recently, Graczyk et al. (2011) 
found 17 EBs in the EMC via the OGEE survey where the eclipse depths changed. In 
two cases the eclipses disappeared altogether. 

The line of apsides of an eccentric binary can rotate due to the effects of General 
Relativity. The ratio of the orbital period of the binary to the rotation period Ggr of the 
line of apsides is given by 


-Pbin 


= 6.35 X 10“^ 


/ Ml +M2 \ 

\a(l -e^)] 


(5) 


where M\ and M 2 are the masses of each star in solar masses, and a is the semimajor 
axis of the relative orbit in solar radii (Kopal 1978). In addition, tidal distortions can 
cause the line of apsides of an eccentric binary to rotate with a period //tide given by 

77 — - 021^21 -I- 022^22 
^tide 
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where ki 2 and k 22 are the so-called internal structure constants of the stars, and where 
the coefficients C 21 and C 22 are given by 


C21 = 

C22 = 



(V) 


In the above equations, and {a> 2 l^K) are the ratios between the actual angular 

rotational velocity of the stars and the rotational velocity corresponding to synchroniza¬ 
tion with the average orbital velocity. The functions f{e) and g(e) are given by 
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(Claret & Gimenez 1993). Note that the expected apsidal period depends on the fifth 
power of the fractional radii of the stars, and that for EBs with periods more than about 
15 to 20 days the apsidal period will usually be factors of several million or more times 
the binary period, assuming the EBs contain solar-type main sequence stars. 

Since the phase difference befween fhe primary and secondary eclipses changes 
wifh lime, “apsidal mofion” due fo General Relalivify or fo lides will give rise lo ETVs, 
where fhe 0-C curves of fhe primary and secondary roughly resemble sine curves fhal 
are 180 degrees ouf of phase. Unforlunalely fhere are no simple close-form expressions 
lo model fhe signal seen in fhe 0-C diagram of an eccenlric EB undergoing apsidal 
mofion. Gimenez & Garcia-Pelayo (1983) give a power series expression up fo fhe fiflh 
power of fhe eccenlricily fhal conlains roughly 50 terms. Eacy (1992) presented an 
exact solution based on an iterative solution of the transcendental equations involved. 


2. The Kepler EB Sample 

The Kepler mission (Borucki et al. 2010) observed over 2000 EBs with a duty cycle on 
the order of 90% or greater, and good photometric precision (Prsa et al. 2011, Slawson 
et al. 2011). These data give us the unique opportunity to measure precise ETs and 
detect deviations from a linear ephemeris. Previously, Gies et al. (2012) measured ETs 
for 41 Kepler EBs using data thr'ough Q9 (« 2 year baseline). They found no evidence 
for short-period companions with periods smaller than w 700 days in the sample, but 
did find evidence for long-ferm frends in 14 systems. Rappaporf ef al. (2013) did a 
much more fhorough survey where fhey measured ETs for 2175 Kepler EBs using dafa 
fhrough Q13 (w 3 year baseline). They identified 39 candidafe friple sysfems. Conroy 
ef al. (2014) presented a cafalog of ETs for 1279 close binaries (e.g. fhe overconfacf and 
ellipsoidal EBs wifh periods generally smaller fhan 1 day) in fhe Kepler sample. There 
were 236 sysfems where fhe ETVs were compafible wifh fhe presence of a fhird body. 

We give an updafe on our own program fo measure accurafe ETs for fhe Kepler 
EBs wifh periods longer fhan abouf 1 day. We have developed a suife of algorifhms and 
codes (collecfively called “TEMPUS”) lo aulomalically measure times of mid eclipse 
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for detached Kepler EBs. TEMPUS was developed on and runs in the Matlab envi¬ 
ronment. An earlier version of the software was described in Steffen et al. (2011). 
TEMPUS has evolved since then, so a new overview is presented here. 

TEMPUS is a computational system for measuring ETs from the Kepler SAP light 
curves. Briefly, a model eclipse profile consfructed from the data is used to find indi¬ 
vidual ETs. Because the ETs may not be described by a simple linear ephemeris, an 
iterative method is used. The iteration begins with an initial set of ET values, usually 
computed from a linear ephemeris. The eclipses in the data are found and locally de¬ 
trended using a cubic spline. The detrending algorithm requires that each eclipse event 
has well-defined “shoulders” so that the first and fourth contact points can be identi¬ 
fied. When such shoulders are presenf, we have found the local detrending works well 
in most cases. Once the individual eclipse events are locally detrended, they are folded 
to produce an eclipse profile. A piecewise cubic Hermite spline (PCHS) model is fhen 
fit to this profile. During this process, nearby data are tested for gaps and other problems 
that would compromise the determination of that ET value. Those ETs that exhibit such 
problems are eliminated from further consideration. The PCHS is then used to improve 
the ET estimates, and those, in turn, provide an improved eclipse profile consfrucfed 
from fhe local defrending and folding. This iteration is run three times. At this point 
the ET estimates are significantly improved and the eclipses with compromised data 
have been eliminated. However, there remains a chance that a good eclipse may have 
been eliminated as the ET estimates were being adjusted. Therefore, using the latest 
ET estimates this iterative process is restarted with any missing ET estimates being ap¬ 
proximated by nearby ones. After three additional iterations, the so-called Pipeline has 
produced very good ET estimates and a very good PCHS model of the eclipse profile. 

The final step in the process is to compute an unbinned PCHS model that when 
binned to the Kepler Eong Cadence exposure time (29.4244 minutes) will best represent 
the eclipse profile made from the detrended and folded data. As an initial approxima¬ 
tion, the PCHS resulting from the six iterations is used. It is binned for each point in the 
latest locally-detrended and folded data and used to improve the ET estimates. These 
ET estimates again determined a revised local detrending and folding of the data with 
a best, unbinned PCHS model being computed and, in turn determining a further, but 
very minor, revision in the ET estimates. The last iteration of this process occurs with 
the final resulf being a PCHS model thaf, when binned, optimally describes the eclipse 
profile produced from the locally detrended and folded data, which are determined by 
the latest ET estimates. 

During this process, ETs that are outliers to the individual model fits are eliminated 
from defining further models but remain as values to be estimated. The error estimates 
for each ET value use the standard approach of sliding the model, in time, across the 
ET estimate and noting when the;t'^ curve rises above the \cr threshold. 

Producing a fully automated process to measure ETs for Kepler EBs is a difficult 
task owing to the wide variety of light curve morphologies (e.g. see Prsa et al. 2011 and 
Slawson et al. 2011): Many binaries are spotted, many of them have pulsations, many 
binaries have deep eclipses, and others are severely diluted. In order to have confidence 
in the results, the TEMPUS pipeline produces many diagnostic plots such as plots of the 
individual 0-C and CPOC diagrams, plots of the observed and model eclipse profiles, 
power spectra of the ETVs, and plots of a small subset of the “raw” light curve. These 
plots were inspected by eye and interesting systems were selected. 
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3. Results 


We used all of the available long cadence data through Q16 as input to TEMPUS. Our 
initial sample consisted of the 1322 binaries from the catalog of Slawson et al. (2011) 
classified as either detached or semidetatched, supplemented with some additional EBs 
that were discovered since that publication. As the orbital period gets shorter and 
shorter, reflection effects and tidal distortion tend to make the out-of-eclipse regions in¬ 
creasingly curved. Since our technique requires well-defined eclipse shoulders fo work, 
we placed a cufoff al orbifal periods shorfer fhan 0.9 days. This cufoff removed a few 
hundred binaries from Ihe sample, leaving 1258. These binaries were fed info TEMPUS 
fo measure Ihe ETs, and 1249 syslems had fheir primary ETs measured successfully. 
The cases fhaf failed eifher lacked well-defined shoulders, or had very low signal-lo- 
noise eclipses. Eccenlric binaries can have primary eclipses and no secondary eclipses, 
so fhe sample of binaries wifh measurable secondary eclipses is somewhaf smaller. In 
Ihe end, a folal of 764 syslems had bolh primary and secondary ETs measured success¬ 
fully. The typical uncerlainlies on Ihe ETs range from a few seconds fo several minutes 
for the noisy cases. 

The diagnostic plots from TEMPUS were inspected and obvious spurious cases 
were removed. Particular attention was paid to the plots showing the final eclipse profile 
and final PCHS model. Based on the quality of the fits to the eclipse profiles, 335 EBs 
were removed from the sample, leaving 914 systems. 




Figure 1. Left: Example CPOC diagrams where the rms scatter in the ETVs is 
smaller than 10 sec. The black points show the primary ETVs and the red points 
show the secondary ETVs. Erom top to bottom, the systems are KIC 9913481, KIC 
10156064, KIC 10191056, and KIC 12356914. Right: The distribution of the rms 
of the primary ETVs in minutes. Most of the systems have reasonably small scatter, 
although the tail of the distribution extends to much higher rms values. 
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Finally, for each binary, the 0-C diagrams for the primary eclipses and secondary 
eclipses (if present) were made. In cases where both primary and secondary eclipses 
were present, an iterative procedure was used to produce a CPOC diagram. These 
plots were visually inspected, and interesting cases were selected. We give below an 
overview of these results. 

3.1. EBs With Small ETVs 

One simple statistic to compute is the rms of the primary ETVs. A large rms in the 
ETVs may indicate the presence of a nearby third body, although one needs to verify 
this on a case-by-case basis (for example, star spots can induce spurious ETV signals 
as discussed below). Eikewise, a small rms in the ETVs might be used to place limits 
on the lack of a third body, although, again, one needs to verify each case individually. 
There may be cases where long-term trends are evident in spite of the small rms, or 
cases where the primary and secondary ETV signals diverge in the CPOC. Eigure 1 
shows four examples where the primary and secondary rms values are smaller than 
about 10 seconds. Eigure 1 also shows the distribution of rms values for the primary 
ETVs. The histogram peaks around 10 seconds. A total of 502 systems have rms values 
smaller than 30 seconds and 741 systems have rms values of smaller than one minute. 



410 420 430 440 

Time (BJD - 2,455,000) 



Figure 2. Left: Eclipses of a star spot in Kepler-47. A dark spot on the primary 
star rotates in and out of view, leading to the variable flux in the out-of-eclipse re¬ 
gions as shown at the bottom. This star spot was partially eclipsed by the secondary 
star during five consecutive primary eclipses. When the eclipse is modeled with a 
symmetric function, a hump appears in the residuals. As the spot moves on the pri¬ 
mary, the phase of the hump in the residuals changes. Right: The O-C values of the 
primary eclipses in Kepler-47 vs the local light curve slope. The eclipses are system¬ 
atically late when the out-of-eclipse flux slopes downward during a primary eclipse 
(this occurs when the dark spot is rotating into view), and systematically early when 
the out-of-eclipse flux slopes upward during a primary eclipse (this occurs when the 
dark spot is rotating out of view). Figures from Orosz et al. (2012). 


3.2. EBs with Star Spots 

In an EB with immaculate stars, the eclipse profiles should be smooth and symmetric 
in time for EBs with small eccentricities. If a star spot (either dark or bright) appears 
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Time (BJD-2,455,000) 


Time (BJD-2,455,000) 


Figure 3. Example CPOC diagrams where the ETVs are affected by star spots, 
giving rise to “random” walk signals (the color scheme is the same as the left panel 
in Eigure 1). Erom top to bottom and left to right, the systems are KIC 4908495, 
KIC 6697716, KIC 6706287, KIC 6863840, KIC 7374746, KIC 9005854, KIC 
10346522, and KIC 12418816. 


on the star that is being eclipsed, then the eclipse profile will appear distorted. The 
distortion can be relatively large in cases where the body in front is much smaller than 
the body that is being eclipsed. If one is using a symmetric model profile fo find the 
time of minimum light, there will be a hump in the residuals of the profile fil and fhe 
resulfing time will have a sysfemafic error. The primary eclipses in Kepler-47 illuslrafe 
fhis quife nicely as shown in Figure 2 (faken from Orosz ef al. 2012, see also Sanchis- 
Ojeda el al. 2012). One or more dark spols on fhe primary rofafe info and ouf of view, 
leading fo a modulation in fhe ouf-of-eclipse regions of the light curve. When the light 
curve has a downwai'd slope near primary eclipse, a dark spot is rotating into view, and 
when the light curve has an upward slope near primary eclipse, there is a dark spot 
rotating out of view. Figure 2 shows a sequence of five primary eclipses where fhe local 
slope of fhe lighf curve goes from being positive fo negative during fhe eclipse. The 
spof is in a differenf place on fhe primary sfar during successive primary eclipses, and 
as a resulf fhe phase of fhe hump in fhe residuals changes. For Kepler-47 and many 
ofher EBs, fhere is a correlation befween fhe ETV of an eclipse lime and fhe local slope 
of fhe lighf curve during lhaf eclipse (Sanchis-Ojeda ef al. 2012), as shown in Eigure 2. 
The presence of fhis correlation in Kepler-47 is a good indication that the signal seen in 
the ETVs of the primary star is spurious, as the ETVs, corrected using this correlation, 
show no signal (Orosz et al. 2012). 

Because star spots seem to come and go, an EB with spots may exhibit a “random 
walk” signal in the CPOC diagram. Such a signal could have systematic deviations 
much larger than the nominal uncertainties in the ETs, but will not remain coherent 
over the long term. During the visual inspections of the CPOC diagrams, 117 EBs with 
random walk signals were identified, and Eigure 3 shows eighf example cases. 
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Figure 4. Left: Example CPOC diagrams for systems with obvious periodicities, 
where the color scheme is the same as in Figure 1. From top to bottom, the EBs 
are KIC 6545018 (Pbin = 3.99 d, = 90.4 d), KIC 9714358 (Pbin = 6.47 d. 
Pout = 103.8 d), KIC 9451096 (Pbi„ = 1.25 d. Pout = 106.9 d), and KIC 5095269 
(Pbin = 18.6 d. Pout = 117.9 d). Right: Two EBs with periodic signals in their CPOC 
diagrams where LTT model hts (blue lines) yield very high tertiary masses. From 
top to bottom, the EBs are KIC 4940201 (Pbin = 8.82 d. Pout = 363.7 d, f(M) - 
4.53 Mq), and KIC 5384802 (Pbin = 6.08 d. Pout = 255.2 d, f(M) = 13.92 Mq). 


3.3. EBs with Obvious Periodicities in the CPOC 

During our visual inspections of the 0-C and CPOC diagrams, we identified 55 EBs 
with obvious periodicities in the primary ETVs, including 7 with periods smaller than 
200 days. Eigure 4 shows 0-C and CPOC diagrams for the four shortest-period ones, 
with ETV periods of 90.4, 103.8, 106.9, and 117.9 days. These periods almost certainly 
indicate the period of the third body. A literature search found mention of three EBs 
that were known before the launch of Kepler with third bodies having periods smaller 
than 200 days. 

The EBs where the primary and secondary ETVs are periodic and track each other 
in the CPOC diagram can be modeled with an ETT orbit (Equation 2) and the minimum 
mass (hereafter the “mass function”) of the third body can be computed using Equation 
4. A few EBs have implausibly large mass functions, and Eigure 4 shows the CPOC 
diagrams of the two EBs with the largest mass functions (4.53 Mq and 13.9234©). As 
shown by Borkovits et al. (2003) and Rappaport et al. (2013), dynamical effects can 
sometimes produce ETV signals that can mimic ETVs due to pure ETT effects. Thus 
one should use extreme caution when modeling ETV signals with a simple ETT model. 

3.4. EBs with Large ETVs, Changing Eclipse Depths, and Tertiary Eclipses 

We have found about a dozen systems where the range on the ETVs alone rules out 
simple ETT effects as the sole cause of the variations. Eigure 5 shows the eight sys- 
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Time (BJD-2,455,000) Time (BJD-2,455,000) 

Figure 5. Example CPOC systems where the range in the ETVs is too large to 
be caused by LTT effects only (the color scheme is the same as in Eigure 1). Erom 
top to bottom and left to right, the systems are KIC 5255552 (Aq-c = 13.3 hr), KIC 
7668648 (Aq-c = 8.7 hr), KIC 5653126 (Aq-c = 6.7 hr), KIC 7955301 (Ag-c = 3.6 
hr), KIC 7955301 (Ap-c = 3.6 hr), KIC 5771589 (Aq-c = 3.5 hr), KIC 5003117 
(Ao-c = 2.5 hr), and KIC 8210721 (Aq-c = 2.5 hr). 


terns with the largest ETV spreads. KIC 5255552 has a spread of about 13.2 hours, 
which would require a displacement of the binary by about 100 AU if the ETVs were 
due entirely to ETT effects. We conclude the cause of these ETV signals is largely 
dynamical. 

There are at least 16 EBs that show significant changes in the eclipse depths due 
to dynamical interactions.^ Eigure 6 shows the light curves and CPOC diagrams for 
eight of these depth-changing systems. In all cases, there are significant ETVs, which 
is another indication that the depth changes are due to dynamical interactions. 

Although hundreds of EBs with a tertiary companion were known prior to the 
launch of Kepler, none were observed to have eclipse events due to the third star. KOT 
126 was the first system with tertiary events identified in the Kepler data (Carter et al. 
2011). We have identified 24 EBs thaf show addifional eclipse events due to a third 
body, including the circumbinary planet systems (there are a few other cases where 
there appear to be two unrelated EBs on the same Kepler pixel). Eigure 7 shows eight 
examples. The presence of these additional eclipse events provides very strong con¬ 
straints on the parameters of the system. Modeling these light curves can be difficult, 
but one is then rewarded with extremely precise parameters (e.g. Carter et al. 2011; 
Doyle et al. 2011). 


'There are also many cases where there are spurious changes in the eclipse depths due to changing amounts 
of contamination due to nearhy stars. Different photometric apertures may be used for different Kepler 
Quarters, and these different apertures can result in different amounts of contamination. 
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Time (BJD-2,455,000) 


Figure 6. Example EBs where the depths of the eclipses have changed due to the 
influence of a third body. The normalized Kepler light curves are shown on the left, 
where different colors denote the Kepler Quarter (black for Ql, Q5, Q9, Q13, and 
Q17; red for Q2, Q6, QIO, and Q14; green for Q3, Q7, Ql 1, and Q15; and blue for 
Q4, Q8, Q12 and Q16). The corresponding CPOC diagrams are shown on the right 
with the same color scheme as in Eigure 1. Erom top to bottom, the systems are 
KIC 4769799, KIC 5255552 (this system fell on the bad CCD module during Q5, 
Q9, Q13, and Q17), KIC 5653126 (note the appearance of the secondary eclipses 
starting in Q8), KIC 5731312, KIC 7289157 (note the additional events due to the 
third body), KIC 7668648, KIC 7670617, and KIC 10319590. 
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Time (BJD - 2,455,000) 


Time (BJD - 2,455,000) 


Figure 7. Some EBs that show tertiary events due to a third star. The color scheme 
for the light curves is the same as in Figure 6. From left to right and top to bottom, the 
EBs are KIC 2835289, KIC 2856960, KIC 4150611, KIC 5255552, KIC 6543674, 
KIC 7289157, KIC 7668648, and KIC 7670485. 


3.5. EBs with Long-Term Trends or Diverging CPOCs 

We found 140 EBs with long-term trends in the ETVs. These trends presumably rep¬ 
resent a small part of a signal with a periodicity that is much longer than our w 4 year 
baseline. Eigure 8 shows four examples. The ETVs for the primary usually, but not 
always, track the ETVs for the secondary in the CPOC diagram. 

We also found 107 EBs where the primary ETVs and secondary ETVs have roughly 
linear trends with opposite slopes in the CPOC diagram. The ETVs with opposite 
slopes occur when only a small part of an apsidal period is observed. Unlike the cases 
with large ETVs or periodic ETVs, having opposite linear trends in the CPOC diagram 
does not necessarily indicate the presence of the third body. As noted earlier, apsidal 
motion can be caused by tidal effects and/or General Relativity. While there are good 
analytic approximations for the rate of the apsidal advance, one needs to know the 
masses to apply the expression for GR precession and the masses and fractional radii to 
apply the expression for the tidal apsidal precession. If one assumes masses near a so¬ 
lar mass for each star', then detailed models of the light curves can give good estimates 
of the eccentricity of the binary and the fractional radii, making the expected rates of 
apsidal motion due to tides and GR computable. We are in the process of modeling 
the light curves to determine what fraction of the EBs where the opposite trends in the 
CPOC diagram could be explained by GR, tidal apsidal motion, or a combination of 
the two. In the absence of detailed light curve modeling, one can make a quick estimate 
by considering only the systems with orbital periods longer than about 20 days, as the 
expected apsidal rates for GR and tides fall off rapidly with increasing period, leaving 
the influence of a third body as the probable cause of the opposite trends in the CPOC 
diagram. 
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Time (BJD-2,455,000) Time (BJD-2,455,000) 


Figure 8. Left: Some EBs that show long-term trends in their ETVs. Erom top to 
bottom, the EBs are KIC 4732015, KIC 5513861, KIC 8429450, and KIC 8553788. 
Right: Some EBs that have rougly linear trends with opposite slopes in the CPOC di¬ 
agram. Erom top to bottom, the EBs are KIC 3247294, KIC 4544587, KIC 5955321, 
and KIC 8553907. 


On a related note, when the primary ETVs and the secondary ETVs show opposite 
trends in the CPOC diagram, this means that one would get two different periods when 
fitting the primary ETs and the secondary ETs separately. We find 184 EBs where fhe 
primary eclipse period differs from fhe secondary eclipse period by more fhan 3cr. Nof 
all of fhese cases have roughly linear signals in fhe CPOC (for example KID 7289157 
in Eigure 6), and fherefore fhe number of sysfems found by fhis fesf exceeds fhe 107 
sysfems found above. There are 136 EBs ouf of fhe 184 wifh orbifal periods 20 days or 
shorter, so in many of fhese cases fhe period differences could be due fo GR and/or tides. 
Eor example, fhe primary and secondary periods of KIC 4344587 differ by w 901 cr (fhe 
mean period is abouf 2.19 days). Gies ef al. (2012) showed fhaf fhis period difference 
could be plausibly explained by fidal apsidal motion, as bofh sfars in fhis shorf-period 
eccenfric binary have relafively large fraclional radii. On fhe ofher hand, we have 48 
sysfems wifh significanf period differences where fhe mean orbifal period is longer fhan 
20 days, and fhe influence of a fhird body will almosf cerfainly be needed fo explain fhe 
period differences in fhese cases. 

3.6. A Brief Note on the Occurrence Rate of Close Triples 

In their limited survey of 41 EBs, Gies et al. (2012) found 14 EBs with long-term 
trends, which represents 34% of the sample. They argued that “this finding is consistent 
with the presence of tertiary companions among a significant fraction of the targets, 
especially if many have orbits measured in decades”. In their more extensive survey, 
Rappaport et al. (2013) found 39 candidate triple systems and after accounting for the 
types of systems to which their search was sensitive, concluded that at least 20% of all 
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close binaries have third body companions. Conroy et al. (2014) found 236 candidate 
triple systems out of 1279 searched, for a rate of about 18%. 

Although this work is ongoing, we can make some preliminary estimates of the 
occurrence rate of triples in our sample. The 0-C diagrams for the primary eclipses 
were visually inspected, and the ETV signal was assigned one of 6 classes: 1 for flat 
signal, 2 for an obviously periodic signal, 3 for a trend that is concave upwards, 4 for 
a trend that is concave downward, 5 for a “wiggle” (see KIC 5731312 in Figure 6 ), 
and 6 for a random walk due to star spots. The ETV signals for categories 2, 3, 4, 
and 5 will most likely be due to third bodies. There were 203 systems with one of 
these four designations out of the 1249 systems measured, which is about 16%. If 
we use the 914 EBs with acceptable profile fits as the sample size, then the fraction 
is 22%. In addition, the CPOC diagrams for all systems with both measured primary 
and secondary eclipses were inspected, and the ones where the primary ETV signal 
diverged from the secondary ETV signal were flagged. There were 112 such cases. In 
total there are 285 unique systems from both lists, which is 22.8% of the sample using 
a sample size of 1249 or 31% using a sample size of 914. As explained earlier, some 
of the EBs with roughly linear trends with opposite slopes in the CPOC diagram may 
not have third bodies (e.g. apsidal precession due to GR or tides may be sufficient to 
explain the ETVs). Thus, the 22.8% or 31% figures are upper limits. Nevertheless, we 
can say that based on our results so far, we can corroborate the results of Rappaport et 
al. (2013) who concluded that the occurrence rate of triples is on the order of 20%. 


3.7. Models of Individual Systems 

We have begun a program to systematically model the light and velocity curves of 
several EBs with large ETVs and/or tertiary eclipse events. The usual binary light 
curve synthesis codes are not adequate as the positions of the stars are not described 
by simple Keplerian orbits. We modified the EEC code (Orosz & Hauschildt 2000) to 
include a dynamical integrator that solves the Newtonian equations of motion to give 
the positions of the stars at any given time. Given the positions of the stars, their radii, 
their relative flux contributions, and their limb darkening properties, model light and 
velocity curves can be computed. 

One of the more striking EBs is KIC 10319590, where the eclipses disappeared 
after Q4 (see Figure 6 ). Spectra were obtained using the echelle spectrograph on the 
Kitt Peak 4m telescope. Fines from the primary and the tertiary star' were detected. 
Our best-fitting to the light and velocity curves is shown in Figure 9. The three masses 
are Mi = 0.95 Mq, M 2 = 0.71 M©, and M 3 = 0.84 M©. The inner and outer periods 
are 21.31 and 457.6 days, respectively. The outer orbit has a moderate eccentricity 
{e = 0.142) and is inclined by 43.1° relative to the binary orbit. This mutual inclination 
is consistent with the predictions of the Kozai Cycle-Tidle Friction model for close 
binary formation where one expects clusters of mutual inclinations near 40 and 140 
degrees (Mazeh & Shaham 1979; Fabrycky & Tremaine 2007). 

KIC 7668648 shows a gradual increase in the eclipse depths and large ETVs (Fig¬ 
ure 6 ). In addition, the third star transits both stars in the EB and is itself occulted 
by the other stars (Figure 7). Spectra were obtained using the echelle spectrograph on 
the Kitt Peak 4m telescope. The spectra are double-lined where both stars in the EB 
are detected. A good photodynamical model was developed, and we find masses of 
Ml = 0.849 M©, M 2 - 0.808 M©, and M 3 = 0.278 M©, and radii of Ri = I.OIOR©, 
/?2 = 0.881/?©, and /?3 = 0.289/?©. The mean periods are 27.094 days for the inner 


14 


Orosz 




Time (BJD - 2,455,000) 


Figure 9. Top: the normalized light curve of KIC 10319590 (points) with the 
best-fitting photodynamical model (red line). The units of the residuals are parts per 
thousand. Bottom: The radial velocities of the primary (black points) and the tertiary 
star (red points) and the best-fitting model curves. 


orbit and 206.4 days for the outer orbit. Figure 10 shows the best-fitting photodynam¬ 
ical model, and one can see that the transits and occultations are well-fit, as are all of 
the primary and secondary eclipses. There is an odd feature in this EB. The eclipse 
event near day -23 (labeled Pi) is deeper than the eclipse event near day -9 (labeled 
S i), so the former was called a “primary” eclipse. The eclipse event labeled P 52 occurs 
51 orbital periods later than Pi and the eclipse event labeled 5 51 likewise occurs 51 
orbital periods later than S However, S 52 is deeper than P 52 , which means the roles 
of the primary and secondary eclipses have reversed! KIC 7668648 is a very important 
system as it has a relatively low-mass star with a good mass and radius determination 
and also from a dynamical point of view, as it is a relatively compact triple. 
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